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ABSTRACT 

A model of the wavelength-integrated scalar irradiance for a vertically 
homogeneous water column is developed. It runs twenty thousand times faster 
than simulations obtained using full Hydrolight code and limits the percentage 
error to less than 3.7%. Both the distribution of incident sky radiance and a wind- 
roughened surface are integrated in the model. Our model removes common 
limitations of earlier models and can be applied to waters with any composition of 
the inherent optical properties. Implementation of this new model, as well as the 
ancillary information required for processing global-scale satellite data, is 
discussed. This new model is fast, accurate, and flexible and therefore provides 
important information of the underwater light field from remote sensing. 

1.0 INTRODUCTION 

One major goal of satellite ocean color is to assess the role of the ocean in the global carbon 
cycle and to examine the factors that affect global climate change (Hooker et al., 1992). 

Although phytoplankton account for less than 1% of the total global plant biomass, they 
contribute to about 40% of global productivity (Falkowski, 1994) and turn over carbon and 
nutrients more rapidly than terrestrial plants (Summerhayes, 1996). To attain the goal of 
analyzing ocean carbon flux, the amount and distribution of ocean phytoplankton must be 
estimated from space. Many efforts to date have focused on relating the radiative signals to 
pigment concentration (Gordon et al . , 1988; O'Reilly et a/., 1998; Cardereta/., 1999). Maps of 
satellite-derived pigment concentration (e.g., biomass) have been widely used to estimate global 
ocean carbon content and productivity (Longhurst, 1995; Platt et al., 1995). However, rates of 
phytoplankton photosynthesis are regulated by many factors including light availability. Hence, 
an accurate estimation of ocean primary production and carbon flux requires a thorough 
description of the underwater light field. 

A common approach for calculating the underwater light field is based on satellite-derived 
estimates of the diffuse attenuation coefficient (e.g. K^A9Q)) (Yeh et al, 1991). Spectral 
information is formulated based on ratios of water-leaving radiance to i^490). Only an average 
value of Kd over the first optical attenuation depth at one spectral waveband can be obtained 
using this approach. Most phytoplankton-based applications, therefore, relate empirical models 
of water content, such as the satellite-derived chlorophyll concentration (Chi), to the spectral 
value of Kd (Sathyendranath, 1 989). 
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Recently, Liu et al. (2001) discussed the sources of error in a Abased model. They 
proposed an alternative fast and accurate model to calculate the underwater scalar irradiance, 

E 0 ( z ) . Their model was based on a large number of simulations compared to results from 
Hydrolight (Mobley and Sundman, 2001) The model considered both the distribution of incident 
sky radiance and a wind-roughened surface and could be applied to Case 1 waters as well as 
gelbstoff-rich Case 2 waters. However, the inherent optical properties (IOPs) were parameterized 
as a function of Chi with a limited range of 0 to 10.0 mg-m . In addition, the Petzold’s average 
particle phase scattering function was assumed and the spectral range limited to 400 - 700 nm. 
These limitations must be removed to apply the model to estimate the underwater light field from 
remote sensing. 

This work extends the model of Liu et al. (2001) to the development of a new model of the 
wavelength-integrated scalar irradiance for a vertically homogeneous water column. This new 
model removes the original model limitations and yet remains both the advantages of speed and 
accuracy. Comprehensive model-to-model comparisons were made to demonstrate that the new 
model could be applied to waters with any composition of IOPs. Implementation of this new 
model, as well as the ancillary information required for processing the global satellite data are 
discussed below. 

2.0 PARAMETERIZING THE INHERENT OPTICAL PROPERTIES 

The earlier model is mainly limited by the parameterization of IOPs (e.g., Chi and the 
wavelength, X ). We started the new model formulation by reexamining the selection of main 
variables for parameterizing IOPs. 

The numerical optical model is simply a model that converts the IOPs to the apparent 
optical properties (AOPs). Its performance relies on how accurate and efficient the IOPs are 
specified. All IOPs can be derived from two principle IOPs: the absorption coefficient a and the 
phase scattering function (Haltrin, 1999), where y is the scattering angle. The key to our 
model is to construct a look-up table (LUT) for quick reference to the vertical profile of the 
average cosine /7(z) . These two principle IOPs should be selected as the main variables in the 
LUT. However, f3(i//) describes the probability of scattering in the y/ direction, which varies 
from a symmetric distribution (molecular scale) to a heavily peak distribution in the forward 
direction (large particle). To consider a large range of the distribution of , it is necessary to 
find an accurate and efficient way to parameterize /3(y/) . 

An analytical phase function for ocean water was proposed by Fournier and Forand (1994), 
parameterized as a function of the backscattering fraction BF, where 

BF’t- ( 1 ) 

b 

A large range in can be obtained by varying the value of BF. For example, RF=0.01 83 
provides a very good fit to the Petzold’s average particle phase scattering function, while BF= 0.5 
yields the pure water phase function. This function has been integrated into Hydrolight as a 
standard method for specifying J3(y/) (Mobley and Sundman, 2001). 



Based on the work of Liu et al. (2001) a LUT with two variables BF (0.0001 — 0.5) and a> 0 
(0.01 - 0.99) was constructed for quick reference to a set of parameters (B 0 ,B^P,B 2 , Q) used by 
the McCormick five-parameter model (McCormick, 1995) 

— = B 0 +B, exp {-Pg) + B 2 exp (-Qg ) . (2) 

pis) 

Each profile of /7(z) is obtained by assuming a black sky a unit light source in Quad; with 
10 quadrants in the plane of the Sim toward the Sim’s direction (Liu et al., 2001). The different 
distribution of sky radiance and the influence of a wind-roughened surface can both be 
considered by following the approach of Liu et al. (2001). A section of the LUT is given in Table 
1 , where the Pearson correlation coefficient PCC manifests how good the McCormick five- 
parameter model fits the vertical profiles of average cosine. Throughout the entire LUT, the 
average value of PCC is 0.999577 and the lowest value is 0.9801 18. 


Table 1 . A section of the model Look-Up Table (LUT). 













3.0 MODEL VALIDATION 


Various models of IOPs were employed to examine whether the new LUT could provide a 
fast and accurate simulation of E 0 ( z ) for any composition of IOPs in a vertically homogeneous 
water column. 

3 . 1 HISTORICAL MODEL OF CASE 1 WATERS 


A set of bio-optical models of IOPs is widely used for Case 1 waters, including the 
absorption coefficient (Prieur and Sathyendranath, 1981; Morel, 1991) 

a(z; A) = a w (A) + 0.06a’ ( A)Chl 065 + F ■ 0.06 • a’ c (440)0*/ 065 exp(-0.014(T - 440)) , (3) 

the scattering coefficient (Gordon etal., 1983; Morel, 1991) 


b(A) = b w (A) + 0.30 




and the normalized phase scattering function f3 (Mobley, 1 994) 
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a w (A) , b w (A) and a* c (A) are given in the papers of Pope and Fry (1997), Smith and Baker 
(1981), and Morel (1988) respectively. F specifies how CDOM (colored dissolved color matter) 
relates to chlorophyll absorption at a reference wavelength of 440 nm. For Case 1 waters, F is 
frequently set to be a value of 0.2 (Prieur and Sathyendranath, 1981; Morel, 1991). 

Two limitations of this set of models are often questioned. J3 p (ys; A) is calculated from the 
analytic Foumier-Forand phase function (Fournier and Forand, 1994) by specifying for the 
backscattering fraction for particle BF P a value of 0.0183 to mimi c the Petzold’s average particle 
phase scattering function (Mobley and Sundman, 2001). A value of 0.01 83 was noted to be too 
high for Chi greater than 2 or 3 (mg-m' 3 ), as discussed, for example by Morel and Gentili (Morel 
and Gentili, 1991). In addition, the actual value of F may vary within a range around 0.2, even 
for the chlorophyll-dominated Case 1 waters (Liu et al., 1999). To extend the applicability of this 
historical set of bio-optical models of IOPs, BF P and F should be set as a free variable rather than 
a constant. 

3.2 HALTRIN’S MODEL OF CASE 1 WATERS 

Another model of IOPs for Case 1 waters was proposed by Haltrin (1999), which contains 
four components, including pure water, large particles (assumed to absorb like Chi and scatter 
like large particles), CDOM (assumed absorbing but non-scatter) and small particles (assumed 
non-absorbing but scatter). This model expresses spectral absorption, spectral scattering, and 
spectral angular scattering coefficients for each component through the concentration of 


chlorophyll (Haltrin, 1999). It is important to note that two types of particles play a role in the 
scattering and the analytic Foumier-Forand phase function (Fournier and Forand, 1 994) is 
employed by specifying the backscattering fraction for large particles BFi a value of 0.00073 and 
small particles BF S a value of 0.03933. 

3 .3 FOUR-COMPONENT MODEL OF CASE 2 WATERS 

A more general model of IOPs for Case 2 waters exists that is similar to Haltrin’s Case 1 
model. However, the amount of each component is not constrained by Chi, the component of 
small particles represents minerals in waters that can both absorb and scatter light, and the 
backscattering fraction for large particles BFi and small particles BF S are not assumed to be a 
constant value. 

The following comparison was made to examine whether our new LUT provides a fast and 
accurate simulation of E 0 (z) for IOPs described by various models. A total of 20 cases were 
compared by randomly specifying values for all parameters, including the solar zenith angle 9 S , 
cloudiness, surface wind speed V winc ?, chlorophyll concentration CM, CDOM ratio F, the 
backscattering fraction for large particles BFi, the mineral concentration M, the mineral type, and 
the backscattering fraction for mineral particles BF S (listed in the legend of Figure 1). Note that a 
set of absorption and scattering spectrum for various types of mineral, including Bukuta’s brown 
earth (M=l), brown earth (M=2), calcareous sand (M=3), yellow clay (M=4) and red clay (M=5), 
is employed in Hydrolight (Mobley and Sundman, 2001). The wavelength-integrated scalar 
irradiance in the PAR range E 0 PAR (z) was computed at 23 different depths ranging from 0 to 50 

m for each case (Figure 1). The percentage error s % is calculated by 

=. 10 ^" - 1 , ( 6 ) 
where the root-mean-square-error in log 10 scale RMSE ]ogW is defined by 


f (‘Og,, log, f 

RMSE vt = ^ : • (7) 

This analysis provides the same emphasis on underestimates and overestimates. 

A very high correlation (r) as well as a large computational speed ratio (CSR) is obtained 
when comparing the results of our model to the results of a Hydrolight run with elastic scattering. 
The percentage error e % and the maximum relative error e mm are both small. Figure 1 
demonstrates that the new model can be applied to waters with various compositions of inherent 
optical properties. 


Model prediction 


(a) IOP's are specified by the historical model of Case 1 waters 



(b) IOP's are specified by the Haltrin's model of Case 1 waters 
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(c) IOP's are specified by the four-component model of Case 2 waters 
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Figure 1. Wavelength-integrated scalar irradiance in the PAR range E 0PAR (z) 







4.0 CONCLUSION 


Implementation of this new model to estimate the underwater scalar irradiance from remote 
sensing requires the ancillary data of surface conditions and the information of IOP’s of water. 
The information of the ambient optical environment and the surface wind speed can be obtained 
by integrating observations from various satellite sensors. Although it is not possible, to date, to 
make a direct estimation of the phase scattering function /3{y/) from space, recent advancements 
(Lee et al., 2001) have enabled a more accurate retrieval of the absorption coefficient a and the 
back scattering coefficient b h . We can estimate the scattering coefficient b from the satellite- 
derived chlorophyll based on the bio-optical model (Eq. 4) and calculate the backscattering 
fraction BF. With these IOP’s and the ancillary data of surface conditions, this model is able to 
provide a new product of the vertical distribution of the scalar irradiance from remote sensing. 
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